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We study a relativistic charged Dirac particle moving in a rotating 
magnetic field. By using a time- dependent unitary transformation, the 
Dirac equation with the time-dependent Hamiltonian can be reduced to a 
Dirac-like equation with a time-independent effective Hamiltonian. Eigen- 
states of the effective Hamiltonian correspond to cyclic solutions of the 
original Dirac equation. The nonadiabatic geometric phase of a cyclic 
solution can be expressed in terms of the expectation value of the compo- 
nent of the total angular momentum along the rotating axis, regardless of 
whether the solution is explictly available. For a slowly rotating magnetic 
field, the eigenvalue problem of the effective Hamiltonian is solved approx- 
imately and the geometric phases are calculated. The same problem for a 
charged or neutral Dirac particle with an anomalous magnetic moment is 
discussed briefly. 
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In quantum mechanics, the Schrodinger equation, even with a timc-indcpcndcnt Hamil- 
tonian, can be solved analytically only in a few cases. With a time-dependent Hamiltonian, 
the problem is more difficult and fewer examples are well studied. One of the well studied 
examples is a nonrelativistic neutral particle with spin and magnetic moment in a rotating 
magnetic field [1-5] . This simple example has received much attention because of its relevance 
to the problem of geometric phases [6-8] , and also because exact solutions are available. In a 
recent work [9] we have studied a nonrelativistic charged particle moving in a rotating mag- 
netic field, with or without a central potential. The case with a central potential can describe 
the valence eletron of an alkaline atom or that of a hydrogen atom under the influence of the 
external magnetic field. The problem can be treated analytically and some exact solutions are 
available. The purpose of this paper is to extend the previous work to a relativistic charged 
Dirac particle (without a central potential). The main interest of the problem is to see how 
the previous results are changed by the relativistic effect. 

Consider a charged particle with spin 1/2, moving in a rotating magnetic field. The motion 
is described by the Dirac equation. This equation can be written in a form similar to the 
Schrodinger one. Thus the ideas of cyclic motion, nonadiabatic geometric phase, etc. for the 
Schrodinger equation [7] also apply here. The Hamiltonian for this system is of course time 
dependent. As in the nonrelativistic case, we use a time-dependent unitary transformation to 
reduce the Dirac equation to a Dirac-like one with a time- independent effective Hamiltonian. 
This is equivalent to transforming the equation into a rotating frame [10] where the magnetic 
field is static. Thus the effective Hamiltonian in that frame is time independent. However, 
as emphasized in Ref. [9], the transformation is merely a mathematical technique, and the 
effective equation in the rotating frame (which is not an inertial system) does not describe a 
real physical problem. The results derived from this equation are not observable in the rotating 
frame. 

As in the nonrelativistic case, it can be shown that eigenstatcs of the effective Hamiltonian 
correspond to cyclic solutions of the original Dirac equation. The nonadiabatic geometric phase 
of a cyclic solution can be expressed in terms of the expectation value of the component of the 
total angular momentum along the rotating axis, regardless of whether the solution is explictly 
available. For a slowly rotating magnetic field, some of the terms in the effective Hamiltonian 
can be treated as small perturbations and the engenvalue problem of the remaining terms can 
be solved exactly. In this approximation, the geometric phases of the cyclic solutions can be 
calculated explicitly. The above described procedure can also be applied to the Dirac-Pauli 
equation [11] for a charged or neutral particle with an anomalous magnetic moment, except 
that the effective Hamiltonian is too complicated and its eigenvalue problem has not been 
solved. 

We begin with the Dirac equation 



where M and q are respectively the mass and electric charge of the particle, is the vector 
potential describing the external electromagnetic fields. The latter is chosen as 



This vector potential produces the magnetic field B(it) which is chosen to be one with a constant 
magnitude B and rotating about some fixed axis at a constant angle 6b and with a constant 




(1) 



Ao = 0, A{t) 




(t) X r. 



(2) 
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frequency uj. The rotating axis is chosen as the z axis of the coordinate system. The magnetic 
field is therefore 

B(i) Sn(i), n(t) = (sin^sCOsa;t,sin^Bsina;i, cos^b). (3) 

We take > without loss of generality. Note that A(t) also generates a time-varying 
electric field. Thus wc arc indeed dealing with a time-varying electromagnetic field. However, 
the electric field does not enter the Dirac equation directly if the particle has no anomalous 
magnetic moment. We write the Dirac equation in the following form: 

indt"^ = E^t)"^, (4a) 

where 

+ Mc^f, (4b) 

where ol = 7^7. To solve the equation we define the orbit angular momentum (in unit of 
h) I — r X p/h, the spin angular momentum s = S/2 where = ie^^^'^^'^^ /2, and the total 
angular momentum j = 1 -|- s; and then make a time-dependent unitary transformation 

*(t) = Wm{t), (5) 

where 



H{t) = cot 



p - -m 

c 



W{t) = exp{-iujtj,), (6) 

and jz is the z-component of the total angular momentum j. This transformation is a general- 
ization of that used in solving the Schrodinger equation for a neutral [2] or charged [9] particle 
with spin in the rotating magnetic field. It is not difficult to show that 

W^{t)H{t)W{t) = H{0). (7) 

Thus we obtain the following equation for 

= //eff^, (8) 

where the effective Hamiltonian 



= H{0) - fuvjz. (9) 

Since Hes is time-independent, Eq. (8) has the formal solution 

$(i) = C/eff(t)$(0), U,s(t) = eM-iHest/n). (10) 

With the obvious relation ^'(O) = $(0), the time-dependent Dirac equation (4) is formally 
solved as 



^{t) = U{t)^{0), U{t) = W{t)U,sit). (11) 

Since U{t) involves no chronological product, this solution is convenient for practical calcula- 
tions. 
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Now we show that eigenstatcs of the effective Hamiltonian correspond to cyclic solutions 
of Eq. (4) . We take the initial condition 

*i(0) = (12) 

where tpi is an eigenstate of the effective Hamiltonian with eigenvalue and calculate \£'i(T) 
where T = 2t:/uj is the period of the rotating magnetic field. Here for convenience we use 
one subscript i to represent all the quantum numbers that is needed to specify an eigenstate. 
Obviously, Ues{t)ipi — e:K.p{—iEit/?i)(pi, vahd for all t, and W{T)ipi — exp(— i27rj2)(/?j. Because 
we can always expand (pi as a linear combination of the eigenstates of jz, we obtain 

*i(T) = exp{-iEiT/h - i7r)*i(0). (13) 

Hence it is indeed a cyclic solution, and the total phase change in a period is 

Si = -EiT/h - TT, mod 271. (14) 

To determine the dynamic phase, we should calculate 

{H{t)), = {^,{t),H{t)^,{t)) = {^^{0),W^H{t)W^,{0)) = {ip^,H{0)ip^). 

Because H{0) — Hes + tujjz, we have 

{H{t)), = Ei + nuj{3z)i. (15) 

Here {jz)i = {jPhiz^d = lit) : j iit)) is the expectation value of jz in the state ^i(t), and 
it is time independent. Note that {H{t))i is also independent of t. Thus the state ^iit) is 
somewhat similar to a stationary state in a system with a time-independent Hamiltonian. The 
dynamic phase is 

A = -h-' r dt {Hit))i = -EiT/h - 27r{jz)i. (16) 
Jo 

Therefore the nonadiabatic geometric phase is 

7i = - A = -TT + '2'K{jz)i, mod 27r, (17) 

and is determined by {jz)i- This is a relativistic generalization of the result for a nonrelativistic 
neutral [2] or charged particle [9] , and has the same form as the corresponding nonrelativistic 
result. It is valid regardless of whether (pi is explicitly available or not, and is convenient for 
approximate calculations if necessary. 

Our next task is to find the eigenvalues and eigenstates of H^s- Since the effective Hamil- 
tonian is somewhat complicated, we have to make some approximation. We assume that uj is 
small such that the term —hujz in ifgff can be treated as a small perturbation. That is, we are 
considering a slowly rotating magnetic field. In the nonrelativistic case [9], the restriction is 
specifically lo <S \q\B/2Mc, and we have argued that this is in fact a rather loose restriction. 
From the following result for the energy levels we would see that the argument also holds in 
the relativistic case. We thus decompose H^q as 
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where 

H^f, = ca-(p- ^Ao) + Mc27°, (19) 



c 

whoes eigenvalue problem will be solved exactly, and 

H'^s = -fi^Jz (20) 

which will be treated as a small perturbation. In Eq. (19) Aq = A(0). Note that H^fi — H{0). 
It is not difficult to show that 

H^f, = eM-^0B3y)Hls exp(z^Bj,), (21) 

where 

Hl^ = ca • (p - U,) + Mc^°, (22) 

= X r and = — -6(0, 0, 1). This is the Hamiltonian of a relativistic charged 
particle in a static uniform magnetic field along the z axis. 

The eigenvalue problem of H^^^ is easy. We write down the eigenvalue equation 

^effC = E\ (23) 

and denote C, = {u, vY where u and v are two-component spinors and the superscript r denotes 
matrix transposition. In terms of the two-component spinors the above equation takes the form 

c(7-{p- -^ A^i u = {E'' + Mc^)v, 

ca-(p- -A,) v^iE""- M^)u. (24) 



c 

One can solve the first for and substitute it into the second. Then an equation in u alone is 
obtained. It can be solved in the cylindrical coordinates (p, 0, z) . The energy eigenvalues are 

= Eln,mm.± = HiMc^f + i^cKf + \q\Bhc[2n, + |m| + 1 - e{q){m + 2m,)]}V2, (25) 

where we use a single subscript i to represent all the quantum numbers and the sign of the 
energy; kz = 2'nnz/d where d is a length in the z direction for box normalization and — 
0, ±1, ±2, rip = 0, 1, 2, ... is a radial quantum number; m = 0, ±1, ±2, . . . and = ±1/2; 
e(g) is a sign function of q. The corresponding eigenfunctions are given by 

u,{p, <!>, z, sz) = 7V,e-"V/2(^^)H4H(^2^2)^^^^^(^^)^ (26) 



where a — ^J\q\B/2hc, the L|^I are Laguerre polynomials [12], XmA^z) eigenstate of 

with eigenvalue m^; and 

Vi{p, (J), z, Sz) = ^o^j^^2 ^ • [P - -^zj Ui{p, (J), Z, Sz). (27) 
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We do not write down the specific form of Vi because it is lengthy and not necessary for the 
subsequent calculations. We just mention that it consists of two terms. One is proportional 
to Mj, and the other involves the factor e'^'^'^''^^^'^Xms+i{-^z) e**^™+^)'^x„i^_i(s2), depending on 
whether = —1/2 or = 1/2, respectively. It should be remarked that neither m nor 
is a good quantum number. A good quantum number related to them is ruj — m -\- rris which 
is the eigenvalue of j^, a conserved quantity. The reason why we need the quantum number 
nis is that for a given rrij there are two kinds of solutions, corresponding to the two values of 
m.. In the above solutions the normalization constants are 



E? 



r(np+ |m| + 1) 



1/2 



(28) 



Note that the solutions in the above forms are not appropriate when Ef = —Mc^, which may 
happen when Uz = 0, Up = 0, m = e(g)|m|, nis = e{q)\ms\. Indeed, one cannot eliminate v 
from Eq. (24) in this case. Rather, one should eliminate u and solve the resulted equation for 
V. The solution reads 

Voomm.-iP, ^, 5,) = 7V^„^^^_e-"''''/2(ap)l-l-^-^Xm.(s.), (29) 



UQOmms-{P, (P, Sz) = ' (p ~ "^2) ^00mm.-(P, 4>, Z, S^). (30) 

We do not write down the specific form for Uoomms- ^ before. The normalization constant in 
the above solution is 

_ ^/2a 

Vr(H + i) 

Thus the equation (23) is completely solved. The solutions of this equation can also be found 
in Ref. [13], but in different forms. The reason why the solution with a specific energy level 
can have different forms is that the energy levels are degenerate. It seems that the solutions 
in our form are more explicit and convenient. 

Now that Eq. (23) is solved, the eigenvalue problem of H^q become trivial. The eigenfunc- 
tions are 

= exp{-i0Bjy)Ci, (32) 

where Q = {ui,Viy, and the corresponding energy eigenvalues are still given by Eq. (25). We 
will use these ip^ as the approximate eigenfunctions of H^s. Of course the explicit functional 
form of is complicated, but this is not necessary in practical calculations. The lowest 
order corrections to the energy eigenvalues are given by the expectation values of H'^q in the 
approximate eigenstates. The corrected energy levels are 

Ei^ E'^ - mjhw cos 9b. (33) 

Now there is no degeneracy in the quantum numbers. 

Our final task is to work out the geometric phase in a period for the ith state, that is, 
a state with the initial condition (12). Since we have only an approximate result </7° for the 
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eigenstate we can calculate the nonadiabatic geometric phase only approximately. The 
result is 



7j — —rrijflB, mod 27r, 



(34) 



where = 27r(l — cos6b) is the solid angle subtended by the trace of the rotating magnetic 
field. For cu <C \q\B/2Mc, this has the same value as the corresponding nonrelativistic result 
[9]. In the above approximation, it can be shown that 



{^i{t),i^i{t)) — mj{sm 9b cosa;^, sin^g sinut, cos 9b)- 



(35) 



Therefore the total angular momentum precesses synchronously with the magnetic field and 
approximately at the same angle 9b with the rotating axis. Then is also (approximately) 
the solid angle subtended by the trace of the total angular momentum. The geometric nature 
of the result (34) is thus obvious. In the nonrelativistic case, both the orbit and spin angular 
momentum precess synchronously with the magnetic field [9]. Here only the total angular 
momentum does. This is the main difference caused by the relativistic effect. 

Now we turn to the problem of a charged or neutral particle with anomalous magnetic 
moment. The Dirac-Pauli equation is 



\ MC ■ 



* = 0, 



(36) 



where /^a is the anomalous magnetic moment, F^i, = d^Au — duA^, and a^^ = i[7^, ■j^]/2. This 
equation differs from Eq. (1) by the last term in the square bracket. When q — it describes 
a neutral particle, otherwise it describes a charged one. With the previous A^, it can be recast 
in the Hamiltonian form: 

ihdt'if ^H{t)^, (37a) 



where 



H{t) = ca 



1 



A{t) 



+ Mc^j° - fi^j'^'E ■ B{t) + in^-f ■ E{t). 



(37b) 



Note that the electric field E = — c ^dtA now enters the equation directly. We make the 
time-dependent unitary transformation (5-6). With some algebra it can be shown that 



W\t)H{t)W{t) = H{0) 



(38) 



still holds in the present case. Thus the equation for $ has the same form as given by Eqs. 
(8-9), and all subsequent discussions until Eq. (17) arc still valid. Unfortunately, H{(}) is too 
complicated even for a neutral particle and we have not been able to obtain any eigenstate of 
it. 

In conclusion we have considered a relativistic charged particle moving in a rotating mag- 
netic field. The Hamiltonian for such a system is time dependent. By making use of a time- 
dependent unitary transformation, the Dirac equation can be reduced to a Dirac-like equation 
with an effective Hamiltonian which is time independent. In this way we obtain a formal 
solution to the original Dirac equation, which determines the time evolution of an arbitrary 
initial state. The time-evolution operator in this formal solution, unlike that for a general 
time-dependent Hamiltonian, involves no chronological product, and thus is convenient for 
practical calculations. Any solution with one of the eigenstates of the effective Hamiltonian 
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as an initial state is a cyclic solution. The nonadiabatic geometric phase in a period for such 
a solution can be expressed in terms of the expectation value of the component of the total 
angular momentum along the rotating axis. This is an exact relation which holds regardless 
of whether the solution is explicitly available, and is convenient for approximate calculations 
whenever necessary. For a slowly rotating magnetic field, the eigenvalue problem of the effective 
Hamiltonian is solved approximately, and the geometric phases are calculated. The difference 
between the rclativistic results and the corresponding nonralativistic ones is discussed. We 
also briefly discussed the same problem for a relativistic particle with an anomalous magnetic 
moment. 
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